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摘 要 :本 研究 提出 了 适 于 模拟 岩石 类 材料 变形 - 拉 裂 过程 的 三 维 可 开裂 拉 格 朗 日 元 方法 ,其 实质 
是 三 维 拉 格 朗 日 元 方法 与 断裂 力学 中 虚拟 裂缝 模型 的 耦合 方法 ,主要 适用 于 立方 体 单 元 。 以 三 维 
拉 格 朗 上 日 元 方法 为 基础 ,引入 虚拟 裂缝 模型 的 具体 过 程 主要 包括 3 个 关键 步骤 :首先 ,通过 节点 周 
围 单 元 的 应 力 获 取 节 点 的 应 力 ,进而 利用 节点 的 最 大 主 应 力 和 岩石 单 轴 抗 拉 强度 进行 节点 分 离 判 
断 ; 然 后 ,选择 与 节点 最 大 主 应 力 重 直 平 面 最 接近 的 单元 边界 作为 裂缝 扩展 方向 ;最 后 ,引入 虚拟 裂 
颖 模型 模拟 裂 颖 的 靖 生 和 扩展 过 程 。 利 用 该 方法 模拟 了 单 轴 拉 伸 岩 样 和 三 点 弯曲 岩 梁 的 交 形 - 拉 
裂 过 程 ,并 初步 分 析 了 上 型 断 黎 能 、 岩 样 高 度 和 单元 尺寸 的 影响 。 同 时 ,通过 将 本 研究 结果 与 单 轴 
拉 伸 岩 样 峰值 应 力 、 临界 位 移 的 理论 解 和 前 人 位 移 控 制 加 载 条 件 下 三 点 弯曲 岩 梁 的 载荷 -位 移 曲线 
的 数值 解 进行 对 比 , 验 证 了 本 研究 方法 的 正确 性 ,这 为 将 来 的 工作 商定 了 良好 的 基础 。 
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A three-dimensional crackable Lagrangian element method for 
modelling tensile cracking processes of rock-like materials 
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Abstract:A three-dimensional crackable continuous-discontinuous method mainly suitable for cube ele- 
ments 1s proposed to simulate deformation and tensile cracking processes of rocks. [Its essence ls a hybrid 
method of the three-dimensional Lagrangian element method and the fictitious crack model in fracture me- 
chanics. Introduction of the fictitious crack model into the three-dimensional Lagrangian element method 
consists of three key steps. First ,the nodal stress is obtained by averaging stresses of elements around the 
node ,and the nodal separation is assessed according to the maximum principal stress of the node and the 
uniaxial tensile strength. Then ,the element boundary closest to the plane perpendicular to the maximum 


principal stress of the node is selected. Finally ,the fictitious fracture model is introduced to simulate the 
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processes of crack initiation and propagation. Deformation and tensile cracking processes of specimens in 


uniaxial tension and three-point bending are simulated by using the present method. The effects of I-type 


fracture energy , rock specimen height and elements size on the direct tension processes are analyzed. 


Meanwhile ,through comparing the present results with the theoretical peak stress and critical displacement 


of the specimen in direct tension and the numerical load-displacement curve of the specimen in three-point 


bending ,the correctness of the proposed method is verified ,laying a good foundation for future work. 


Key words :three-dimensional crackable Lagrangian element method ;fictitious crack model; uniaxial ten- 


sion ;three-point bending; mesh-dependence ; size effect 


在 高 地 应 力作 用 下 隧洞 围 宕 的 板 裂 化 采矿 活 
动 引 起 的 地 下 岩层 离 层 和 地 表 和 裂缝 中 , 拉 裂 现象 是 
十 分 普遍 的 。 因 此 ,研究 岩石 类 材料 的 变形 - 拉 
列 过 程 具 有 重要 的 理论 和 实际 意义 。 

现 有 的 数值 模拟 方法 可 被 划分 为 3 大 类 :连续 
方法 . 非 连续 方法 和 连续 - 非 连 续 方 法 。 连 续 方 法 适 
于 模拟 连续 介质 的 变形 和 破坏 。 非 连续 方法 适 于 模 
拟 非 连续 介质 的 开裂 和 运动 ”。 连 续 方法 的 不 足 
是 功能 受 限 , 仅 能 在 一 定 程 度 上 处 理 非 连续 问题 。 
非 连续 方法 的 不 足 是 应 力 和 应 变 的 精度 低 , 通 稼 需 
要 在 未 开裂 时 就 引入 了 一 些 刚 度 参数 。 为 了 避免 上 
述 两 种 方法 的 各 上 自 不 足 , 连 续 - 非 连续 方法 (耦合 方 
法 ) 应 运 而 生 ”" , 兼 具 上 述 两 种 方法 的 各 自 优势 。 

二 维 拉 格 明日 元 方法 与 离散 元 方法 的 耦合 方法 
是 一 种 国内 自主 研发 的 连续 - 非 连续 方法 “” 。 该 
方法 适 于 模拟 连续 介质 向 非 连 续 介质 转化 和 非 连 续 
介质 进一步 演化 ,其 正确 性 已 从 多 方面 进行 了 验 
证 。 文 献 L[14] 模 拟 了 单 轴 拉 伸 岩 样 的 变形 - 开 
裂 过 程 ,并 通过 将 分 离 节 点 的 最 大 法 回力 等 的 数值 
解 与 理论 解 对 比 , 验证 了 该 方法 的 正确 性 ; 文献 
[15 模拟 了 三 点 弯曲 崖 架 的 变形 -开裂 过 程 ,并 通过 
将 梁 的 载 傈 -位移 曲线 等 的 数值 解 与 前 人 数值 解 对 
比 ,验证 了 该 方法 的 正确 性 ;文献 L16 | 在 原 有 二 维 
连续 - 非 连 续 方 法 的 基础 上 引入 了 开 型 断裂 能 ,并 模 
拟 了 下 接 剪 切 岩 样 的 变形 -开裂 过 程 ,通过 将 岩 样 上 
辣 面 平均 剪 应 力 -水 平 位 移 曲 线 和 水 平 坐 标 轴 所 围 
面积 与 下 型 断裂 能 取 值 对 比 , 验 证 了 该 方法 的 正确 
性 。 应 当 指 出 ,该 方法 仅 适 于 二 维 问题 ,并 没有 广泛 
的 适用 性 。 因 此 ,迫切 需要 发 展 三 维 连续 - 非 连续 
ts 

本 人 研究 在 三 维 拉 格 度 日 元 方法 的 基础 上 引入 了 
断裂 力学 中 的 虚拟 裂缝 模型 ,提出 了 适 于 模拟 贿 石 
类 材料 变形 - 拉 裂 过 程 的 三 维 可 开裂 拉 格 朗 日 元 方 
法 ,并 模拟 了 单 轴 拉 伸 岩 样 的 变形 - 拉 裂 过 程 , 初 步 
分 析 了 工 型 断裂 能 、 岩 样 高 度 和 单元 太吉 的 影响 。 


还 模拟 了 三 点 弯曲 岩 梁 的 变形 - 拉 裂 过 程 。 本 研究 
工作 为 将 来 工作 (引入 三 维 块 体 接触 算法 ) 奠定 了 
良好 的 基础 。 


1 三 维 可 开 复 拉 格 明日 元 方法 简介 


1.1 三 维 拉 格 朗 日 元 方法 的 基本 原理 


拉 格 明日 元 方法 的 实质 是 基于 拉 格 明日 网 格 的 
连续 介质 有 限 差分 单元 离散 的 中 心 差分 求解 方 
法 ”。 首 先 ,利用 几何 方程 和 本 构 方程 求解 应 力 和 
应 变 ;然后 ,利用 虚 功 原理 求解 节点 的 弹性 力 ;最 后 ， 
利用 中 心 差分 形式 的 运动 方程 求解 节点 的 加 速度 。 
重复 上 述 3 个 步 又 ,直到 获得 满意 的 解 为 止 。 

可 将 三 维 情形 下 连续 介质 离散 成 具有 8 个 市 点 
的 六 面体 单元 。 这 样 , 对 于 位 于 连续 介质 表面 的 1 
个 节点 ,可 能 被 4 个 六 面体 单元 共用 (图 1a 中 市 点 
4 被 单元 ~ 9 共用 ) ,而 对 于 位 于 连续 介质 内 部 的 
1 个 节点 ,可 能 被 8 个 六 面体 单元 共用 (图 1b 中 市 
点 B 被 单元 (D ~ (B® 共用 )。 应 当 指 出 ,图 1 中心 ~(® 
均 为 单元 编号, 分别 被 标记 在 每 个 单元 面 咎 读者 的 
面 的 左上 角 。 有 时 ,为 避免 单元 编号 和 其 他 标记 发 
生 重 合 ,略微 调整 位 置 。 后 文 同 样 存 在 类 似 情况 ,不 


(a) 位 于 连续 介质 表面 的 节点 4 
图 1 位 于 连续 介质 表面 和 内 部 的 节点 


Fig.1 Nodes located on the surface and inside 


(b) 位 村 连续 介质 内 部 的 节点 刀 


of the continuous medium 
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1.2 三 维 可 开裂 拉 格 朗 日 元 方法 


三 维 可 开裂 拉 格 明日 元 方法 的 实质 是 三 维 拉 格 
明日 元 方法 与 断裂 力学 中 的 虚拟 裂缝 模型 的 看 合 方 
法 ,其 3 个 关键 步骤 如 下 。 

1) 市 点 分 离 判 断 。 在 三 维 拉 格 朗 日 元 方法 中 ， 
节点 本 身 并 无 应 力 , 只 有 单元 才 有 应 力 。 目 然 地 , 符 
要 判断 某 一 厄 操 是否 发 生 分 离 , 应 考虑 该 节点 周转 
单元 的 应 力 。 节 点 7 的 应 力 0; 为 该 节点 周围 单元 的 
应 力 的 平均 值 , 即 


Fe 
= Lo (1) 
m=1 


式 中 :n 为 共用 该 节点 的 单元 数目 ;cy” 为 nn 个 单元 
中 任 一 个 的 全 应 力 ,m =1,2,…,n。 

在 三 维 拉 格 朗 日 元 方法 中 ,采用 了 混合 离散 方 
法 以 避免 沙漏 。 先 采用 两 种 方式 对 同一 单元 进行 离 
散 , 即 两 种 覆盖 (图 2 中 1 ~8 为 六 面体 单元 各 个 节 
点 的 编号 ) ,再 对 两 种 结果 进行 平均 。 这 样 ,o,” 为 

: 之 of VA® 
人 (2) 
> TAR 

式 中 :V'" 为 某 一 六 面体 单元 被 混合 离散 成 的 10 个 
四 面体 子 单元 中 任 一 个 的 体积 ;ec 六 为 任 一 四 面体 
子 单元 的 应 力 。 


图 2 两 种 覆盖 
Fig.2 Two kinds of overlays 


利用 最 大 拉 应 力 准 则 ”对 节点 是 否 发 生 分 离 
进行 判断 , 即 


O03 > Or) (3) 

式 中 :o, 为 岩石 的 单 轴 抗 拉 强 度 ;o3 为 最 大 主 应 力 。 
2 ) 分 离 方 向 选择 。 市 点 的 分 离 方 同 与 岩石 的 开 
裂 路 径 紧 密 相 关 。 目 前 , 仅 允 许 沿 单元 边界 开裂 。 
显然 ,模拟 复杂 问题 时 ,单元 尺寸 不 应 过 大 。 否 则 ， 
可 能 导致 饥 齿 形 开 裂 路 径 。 应 当 指 出 , 裂 颖 的 理论 
开裂 方向 可 能 会 位 于 单元 内 部 ,但 在 本 人 研究 方法 中 
裂缝 并 不 能 沿 着 单元 内 部 扩展 ,而 是 选择 与 理论 开 
裂 方 癌 最 接近 的 单元 边界 作为 实际 开裂 方向 ,这 具 
有 无 需 改 变 时 步 长 度 的 特点 。 对 于 以 下 研究 的 正六 
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面体 单元 和 加 载 条 件 , 不 存在 锯齿 形 开裂 路 径 的 可 
能 性 ,因为 岩石 的 拉 腹 发 生 在 水 平 或 垂直 方向 上 。 

选择 与 节点 的 os 垂直 的 平面 M' (该 平面 过 单 
元 (D ~(@ 共 用 的 节点 1 , 且 垂 直 于 节点 1 的 os , 见 
图 3a, 图 3b 为 图 3a 的 主 视 图 ) 最 接近 的 单元 边界 作 
为 裂缝 扩展 方向 。 在 此 分 离 方 向 上 ,一 个 市 点 将 分 
离 成 两 个 分 离 节 点 ,其 中 一 个 分 离 节 点 继承 原 节 后 
编号 ,而 另 一 个 分 离 记 点 拥有 新 编号 。 例 如 ,在 节点 
发 生 分 离 前 ,单元 山 和 字 共 用 节点 1 和 2( 图 4a) ;在 
节点 发 生 分 离 后 ,单元 中 的 节点 和 2 为 继承 原 方 
点 编号 的 分 离 节 点 ,单元 巴 的 节点 3 和 4 为 拥有 新 
编号 的 分 离 节 点 ( 见 图 4b) 。 


| ABCOD) 
AI 

区 4 ©(®) 
0O(2) X 


(a) 与 节点 0 垂直 的 平面 MI (b) XYO7 平 面 的 主 视 网 
图 3 分 离 方 向 确定 示意 图 


Fig.3 Schematic of determination of the separating direction 


(a) 点 发 生 分 离 前 


(b) 节点 发 生 分 离 后 
图 4 市 点 分 离 过 程 
Fig.4 The separating processes of nodes 
3) 虚拟 裂 缝 模型 引入 。 在 常规 的 有 限 元 方法 或 
二 维 拉 格 明日 元 方法 中 引入 虚拟 裂缝 模型 的 成 功 案 
例 已 有 不 少 “” 。 与 常规 的 有 限 元 方法 相 比 ,本 研 
究 方 法 适 于 处 理 大 变形 和 塑性 流动 问题 。 在 本 人 研究 
中 ,虚拟 裂缝 模型 的 引入 是 通过 在 分 离 节 后 上 施加 
黏 聚 力 引 起 的 节点 力 来 实现 的 。 裂 颖 的 张 开 必 然 要 
消耗 能 量 ,单位 面积 上 消耗 的 能 量 即 为 断裂 能 。 由 
于 荆 型 断裂 能 一 般 不 容易 测量 且 岩 石 宏观 藤 裂 的 细 
观 机 制 之 一 是 拉 裂 ,所 以 只 引入 了 工 型 断裂 能 G1。 
随 着 裂缝 的 逐渐 张 开 , 裂 颖 承载 能 力 必然 逐渐 减弱 ， 
即 阻碍 询 缝 张 开 的 力 ( 询 颖 面 颖 凤 力 引起 的 市 扩 
力 ) 必然 减少 。 在 Gt 未 被 完全 消耗 之 前 , 裂 颖 仍 具 
有 一 定 的 承载 力 , 此 时 的 裂 颖 为 虚拟 裂 颖 。 否 则 , 裂 
颖 为 真实 裂 颖 。 在 理论 上 ,对 两 个 分 离 节 点 施加 的 
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黏 聚 力 引起 的 节点 力 应 大 小 相等 ,方向 相反 ,因为 它 
们 是 一 对 内 力 。 黏 聚 力 引 起 的 节点 力 ,与 裂缝 张 
开 度 w 的 关系 为 
WL ,0<w<w, a 
0 ， WwW >w 
式 中 :w ,为 w 的 最 大 值 ,w, =2G1/o,; 下 为 节点 发 生 
分 离 时 分 离 节点 的 弹性 力 。 
应 当 指 出 ,由 于 只 对 岩石 拉 裂 问题 进行 模拟 ， 
此 上 述 竺 聚 力 不 包 括 切身 黏 聚 力 。 另 外 , 黏 聚 力 乘 
以 作用 面积 即 为 黏 聚 力 引 起 的 节点 力 。 


2 算 例 及 结果 分 析 


2.1 单 轴 拉 伸 岩 样 


2.1.1 计算 模型 与 方案 

岩 样 被 划分 成 略 干 正六 面体 单元 ,在 岩 样 上 端 
面 施 加 垂直 向 上 的 速度 v, 其 大 小 为 0.01 m/s, 在 宕 
样 下 端面 施加 深 动 支 座 , 设 岩 样 高 度 为 ,在 H/2 
处 , 即 岩 样 上 下 对 称 面 上 设置 允许 开裂 位 置 ( 图 5)。 
计算 在 大 变形 条 件 下 进行 ,不 计 重力 。 部 分 计算 参 
数 如 下 :密度 为 2700 kg/m ,弹性 模 量 为 20 GPa, 泊 
松 比 为 0.2, oo, 为 2 MPa。 

本 人 研究 共 设 计 7 个 方案 。 方 案 1 ~ 方案 3 的 差 
别 仅 在 于 Gr 的 取 值 , Gt 分 别 为 100 .300 .500 N/m， 
岩 样 尺寸 均 为 0.5 m x0.5 m x0.5 m, 单 元 尺寸 均 为 
0.0625 m。 方案 1、 方案 4 和 方案 5 的 差别 仅 在 于 
单元 尺寸 ,后 两 者 的 单元 尺寸 分 别 为 0. 125 m 和 
0.031 25 m。 方案 1 方案 6 和 方案 7 的 差别 仅 在 于 
宕 样 高 度 ,后 两 者 的 岩 样 高 度 分 别 为 0.25 m 和 
0.73 m, 


允许 开裂 位 轩 


图 5 计算 模型 
Fig.S Calculation model 
2.1.2 结果 分 析 
图 6 给 出 了 方案 1 ~ 方案 3 的 载 集 - 位 移 曲线 。 
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其 中 ,此 载 谷 为 岩 样 上 疾 面 的 应 力 乘 以 原始 横 和 截面 
面积 ,此 位 移 为 岩 样 上 端面 任 一 节点 的 垂直 位 移 。 
由 图 6 可 以 发 现 以 下 结论 。 


载荷 /10 N 


位 移 /10”m 
图 6 方案 1 ~ 方案 3 的 载 符 - 位 移 曲 线 
Fig.6 Load-displacement curves of schemes 1 ~3 

1) Gt 对 岩 样 的 载荷 -位 移 曲线 的 峰 后 斜率 有 明 
显 的 影响 。 随 着 6 的 增加 , 峰 后 斜率 变 大 , 即 载荷 - 
位 移 曲 线 的 峰 后 部 分 变 平缓 。 

2) Gt 对 岩 样 的 峰值 载荷 影响 甚 微 。 方 案 1 ~ 
方案 3 的 峰值 载荷 分 别 为 5.085 x10°、5.086 x10”、 
5.088 x 105 N, 对 应 的 峰值 应 力 分 别 为 2.034、 
2.034 4 、2.035 MPa, 这 与 宕 样 的 理论 抗 拉 强度 , 即 计 
算 中 所 采用 的 抗 拉 强 度 ( o，=2 MPa) 极 为 接近 。 这 
在 一 定 程度 上 说 明了 当前 计算 结果 的 正确 性 。 应 当 
指出 ,3 个 方案 的 峰值 衔 载 略 有 差异 ,这 是 因为 显 式 
差分 方法 造成 的 误差 。 

3) 6 对 峰 后 载荷 刚 降 至 零 时 的 位 移 ( 临界 位 
移 ) 有 明显 影响 。 随 着 6# 的 增加 ,临界 位 移 显 著 增 
加 。 方案 1 ~ 方案 3 的 临界 位 移 分 别 为 1. 01 x 
10“3x10“、4.9 x10“m。 应 当 指 出 ,该 临界 位 
移 与 虚拟 裂缝 的 最 大 法 癌 张 开 度 相对 应 。 此 后 , 虚 
拟 裂 颖 将 成 为 真实 裂 颖 (不 再 具有 承受 拉 伸 载 傈 的 
能 力 ) 。 当 虚拟 裂 颖 处 于 向 真实 裂 颖 转化 的 临界 状 
态 时 , 兰 样 中 各 单元 将 处 于 零 应 力 状 态 (无 应 力 和 变 
形 ) 。 此 时 , 岩 样 上 端面 任 一 节点 的 位 移 将 完全 由 真 
实 裂 颖 所 贡献 。 这 也 在 一 定 程 度 上 说 明了 当前 结果 
的 正确 性 。 应 当 指 出 ,在 满足 数值 稳定 的 前 提 下 ,为 
了 提高 计算 效率 ,时 间 步 长 取 值 略 大 ,这 使 岩 样 允许 
开裂 位 置 的 所 有 市 点 同时 发 生 分 离 。 

图 7 ~8 分 别 给 出 了 方案 1 和 方案 3 的 变形 - 拉 
裂 过 程 。 其 中 ,单元 颜色 代表 co ,大 于 零 代 表 拉 伸 ， 
本 人 研究 中 其 他 云图 也 是 如 此 ;市 点 位 移 的 放大 倍数 
为 10, 这 是 为 了 方便 显示 真实 裂缝 ;NN 为 时 步 数 目 ， 
其 与 岩 样 上 瑞 面 任 一 市 点 的 慌 直 位 移 或 加 载 时 间 成 
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正比 。 由 此 可 以 发 现 : 随 着 NN 的 增加 , 宕 样 中 各 单元 ”被 拉 裂 成 上 下 两 部 分 ,此 后 , 兰 样 中 各 单元 的 os 随 着 
的 os 先 增加 (将 样 的 变形 处 于 峰 前 阶段 ) ,直至 宕 样 ”NN 的 增加 而 下 降 ( 岩 样 的 变形 处 于 应 变 软化 阶段 ) 。 


x10° 
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图 7 方案 1 的 变形 - 拉 裂 过 程 中 os 的 时 空 分 布 


Fig.7 Spatiotemporal distributions of os during the deformation and tensile cracking process in scheme 1 
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图 8 方案 3 的 变形 - 拉 裂 过 程 中 os 的 时 空 分 布 

Fig.8 Spatiotemporal distributions of og; during the deformation and tensile cracking process in scheme 3 

图 9 给 出 了 方案 1 方案 4 和 方案 5 的 载 傈 -位 

移 曲 线 ,图 10 给 出 了 方案 5 的 变形 - 拉 裂 过 程 。 可 
以 发 现 , 载 和 荷 -位 移 曲 线 的 峰 后 部 分 的 网 格 依赖 性 几 
图 11 给 出 了 方案 1 方案 6 和 方案 7 的 载 傈 -位 

移 曲 线 ,图 12 给 出 了 方案 7 的 变形 - 拉 裂 过 程 。 由 
此 可 以 发 现 , 载 答 - 位 移 曲线 的 峰 后 斜率 明显 受到 岩 
样 蜗 度 的 影响 。 宕 样 高 度 越 大 ,和 载 集 -位 移 曲线 的 峰 


载荷 /10 N 


后 部 分 越 陡峭 。 也 就 是 说 ,载荷 -位 移 曲 线 的 峰 后 部 ee 
分 明显 具有 尺寸 效应 ,其 斜率 不 能 被 视 为 材料 参数 图 9 方案 1 方案 4 和 方案 5 的 载荷 -位 移 曲线 
(本 构 参数 ) ,而 是 一 种 由 虚拟 裂 颖 和 之 外 弹性 体 构 Fig.9 Load-displacement curves of scheme 1 ， 
成 的 结构 的 啊 应 。 scheme 4 and scheme 5 
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图 10 方案 5 的 变形 - 拉 裂 过 程 中 os 的 时 空 分 布 


Fig. 10 Spatiotemporal distributions of og; during the deformation and tensile cracking process in scheme 5 
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位 移 /10”m 


图 11 方案 1 方案 6 和 方案 7 的 载荷 -位 移 曲 线 
Fig. 11 


Load-displacement curves of scheme 1 ， 


scheme 6 and scheme 7 
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(ec) N=4 500 


(d) N=5 100 
单位 ;Pa 
图 12 方案 7 的 变形 - 拉 裂 过 程 中 o 的 时 空 分 布 
Fig. 12 Spatiotemporal distributions of og; during the 
deformation and tensile cracking process in scheme 7 
另外 ,上 述 3 个 方案 中 均 只 有 一 条 虚拟 裂缝 ,而 
且 临 界 位 移 和 峰值 载荷 基本 相同 ,这 说 明了 3 个 方 
案 消耗 的 能 量 ( 载 和 荷 -位 移 曲 线 与 横 轴 所 围 成 的 面 
积 ) 基 本 相同 。 单 位 面积 上 所 消耗 的 能 量 即 为 G1， 
3 个 方案 的 Gt 完全 相同 ,这 显然 是 正确 的 。 
2.2 三 点 弯曲 岩 梁 
2.2.1 计算 模型 与 方案 
岩 染 尺寸 为 10 mm x3 mm xl mm ,被 划分 成 40 x 
12 x4 个 正六 面体 单元 。 单 元 尺寸 为 2.5 x10“ m。 
如 图 13 所 示 , 岩 架 最 初 水 平 放置 ,内 染 焉 度 方 向 为 
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2 方 癌 , 岩 染 高 度 方 加 为 了 方向 。 在 知 染 的 左右 侧 
面 的 下 边 分 别 施 加 球 贸 链 和 滚动 文 座 。 在 岩 染 的 上 
端面 且 位 于 路 中 位 置 的 节点 上 施加 垂直 辐 下 的 其 
大 小 为 0.01 m/s。 设 宕 梁 的 跨度 为 上 ,在 宕 染 的 跨 
中 位 置 设置 允许 开裂 位 置 (被 施加 速度 的 节操 不 允 
许 开裂 )。 计 算 在 大 变形 条 件 下 进行 ,不 计 重力 。 部 
分 计算 参数 如 下 :密度 为 2700 kg/m ,弹性 模 量 为 
100 MPa, 泊 松 比 为 0.167,o, 为 1 MPa。 


图 13 ”计算 模型 
Fig. 13 Calculation model 

本 人 研究 共 设 计 3 个 方案 。 方 案 1 ~3 的 差别 仅 
在 于 G1 的 取 值 ,分 别 为 100 .200 .300 N/m。 
2.2.2 结果 分 析 

图 14 给 出 了 方案 1 ~3 的 载荷 -位 移 曲 线 , 同 
时 ,还 给 出 了 文献 [20|] 的 数值 解 。 其 中 ,此 载 集 为 
岩 梁 跨 中 位 置 被 施加 速度 的 节点 的 弹性 力 , 此 位 移 
为 上 述 任 一 节点 的 垂直 位 移 。 


位 移 /10”m 
图 14 方案 1 ~3 和 文献 [20] 的 载荷 -位 移 曲 线 
Fig. 14 Load-displacement curves of schemes 
1 ~3 and reference [20] 

图 15 ~ 图 16 分 别 给 出 了 方案 1 和 方案 3 的 变 
形 - 拉 裂 过 程 。WN 与 岩 梁 跨 中 位 置 被 施 加 速度 的 市 
点 的 垂直 位 移 或 加 载 时 间 成 正比 。 巾 图 15 和 图 16 
可 以 发 现 ,o3; 总 集中 在 裂缝 尖端 附近 , 随 着 裂缝 的 
不 断 间 上 扩展 ,os 集中 区 不 断 问 上 迁移 。 

由 图 14 ~ 图 16 可 以 得 出 以 下 结论 。 
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过 
1) Gr 对 岩 梁 的 峰值 载荷 和 载荷 -位 移 曲 线 的 峰 
后 和 斜 变 有 明显 的 影响 。 随 着 6# 的 增加 ,载荷 峰值 的 
绝对 值 增加 ,例如 ,方案 1~3 的 载 衙 峰值 分 别 为 
-0.9564、-0.9985、-1.007N; 载 荷 -位 移 曲线 的 
峰 后 斜率 变 大 。 这 与 文献 121 | 的 结果 在 定性 上 
相符 。 
2) 在 峰值 载荷 之 前 ,三 点 弯曲 岩 梁 已 发 生 了 拉 


(a) N=6 000 


(d) N=200 000 


(b) N=130 000 


(e) N=220 000 


裂 。 在 不 同方 案 中 ,峰值 载 傈 时 的 拉 裂 程度 有 所 不 
同 ,例如 ,方案 1 ~3 的 拉 裂 颖 长 度 分 别 为 0.001 75、 
0.002 25 .0.002 3 m。 

3) 文 献 [20] 利 用 扩展 有 限 元 法 模拟 了 同一 问 
题 ,方案 1 和 文献 [20 |] 的 载 答 -位 移 曲 线 基本 吻合 。 
应 当 指 出 ,方案 1 的 宕 梁 尺 寸 , 弹 性 模 量 、 泊 松 比 和 
Gr 等 计算 参数 均 与 文献 [20] 的 一 致 。 
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图 15 方案 1 的 变形 - 拉 裂 过 程 中 os 的 时 空 分 布 


Fig. 13 Spatiotemporal distributions of os during the deformation and tensile cracking process in scheme 
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图 16 方案 3 的 变形 - 拉 裂 过 程 中 os 的 时 空 分 布 


Fig. 16 Spatiotemporal distributions of gs during the deformation and tensile cracking process in scheme 3 
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3 结 论 


1) 本 人 研究 提出 了 适 于 模拟 岩石 类 材料 变形 - 拉 
裂 过 程 的 三 维 可 开裂 拉 格 明日 元 方法 ,其 实质 是 三 
维 拉 格 明日 元 方法 与 断裂 力学 中 虚拟 裂 颖 模型 的 耦 
合 方法 ,主要 适用 于 立方 体 单元 。 以 三 维 拉 格 半日 
元 方法 为 基础 ,引入 虚拟 裂缝 模型 的 具体 过 程 主要 
包括 3 个 关键 步骤 :首先 ,进行 节点 分 离 判 断 ; 然 后 ， 
确定 裂缝 扩展 方向 ;最 后 ,引入 虚拟 裂 颖 模型 。 

2) 对 于 单 轴 拉 伸 儿 样 的 载 谷 -位 移 曲 线 , 随 着 I 
型 断裂 能 的 增加 ,临界 位 移 增 加 , 峰 后 斜率 变 大 , 峰 
值 载 集 不 受 影响 ; 峰 后 部 分 具有 尺寸 效应 且 几 乎 不 
具有 网 格 依赖 性 。 对 于 三 点 弯曲 岩 深 的 载荷 -位 移 
曲线 , 随 春 工 型 断裂 能 的 增加 , 峰 后 斜率 变 大 ,峰值 
载体 的 绝对 值 增加 。 同 时 ,通过 将 本 人 研究 结 末 与 单 
轴 拉 伸 知 样 峰值 应 力 、 临 界 位 移 的 理论 解 和 前 人 位 
移 控 制 加 载 条 件 下 三 点 斌 曲 宕 梁 的 载 从 -位移 曲线 
的 数值 解 进 行 对 比 , 验 证 了 本 人 研究 方法 的 正确 性 ,这 
为 将 来 的 工作 问 定 了 良好 的 基础 。 
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